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LONGITUDINAL MOTIONS AND LOADS AS OBTAINED FROM 
FLIGHT MEASUREMENTS OF A SWEPT-WING BOMBER 
By James J. Donegan and Carl R. Huss 


SUMMARY 


Some measurements obtained during elevator-pulse maneuvers with a 
Boeing B-^7A airplane in the Mach nimiber range from 0.62 to 0.82 are 
analyzed to obtain frequency responses in pitching velocity and normal 
load factor up to frequencies of ^5 radians per second. Data -reduct ion 
methods are applied to these measurements to yield transfer coefficients 
and longitudinal stability derivatives. 

The effects of stinictural flexibility are examined by considering 
the B-^7A airplane configuration with three different amounts of flexi- 
bility and comparing the maximum responses of each to the same inputs. 

The ''elastic” airplane contains all the flexibility characteristics of 
the B-47A as obtained from flight data^ the "quasi-steady" airplane has 
the frequency response of the elastic airplane with the structural modes 
removed^ and the "rigid" airplane has a frequency response with all 
effects of structural flexibility removed. Load factor and pitching 
velocities beyond the test range are calculated by mathematical extensions 
of the data to more hazardous conditions. Maximimi time responses for the 
various airplanes subjected to triangular elevator pulses are then pro- 
portioned with one another to obtain quantitative measures of the effects 
of flexibility. The effects of including or omitting phasing between the 
various modes are also studied. 


INTRODUCTION 


In many investigations of the loads on relatively rigid airplanes^ 
it is possible by means of a few well-chosen flight tests in a given speed 
range to obtain data which can be used to predict the loads which would 


2 


KA.CA RM L3kU.6 


occur under more hazardous conditions. For instance, hy proper analysis 
of the flight data, effective airplane parameters can he determined and 
the loads following much more severe control motions at any altitude can 
he predicted. References 1 and 2 are among several reports Tihich present 
methods of reducing data to the required coefficient form. 

Although the effects of small amounts of flexibility are handled hy 
using the above methods and terming the reduced parameters as effective 
values, the effects of large amounts of flexibility, especially in com- 
bination with swept-wing configurations, cannot be bypassed in this 
manner. In such cases, structural modes and structural deformations which 
cause aerodynamic-center shifts and changes in aircraft frequency and 
damping became important variables. Usxially these effects of additional 
structural modes on aircraft loads can be reasoned qualitatively but can- 
not be predicted numerically. 

The investigations on loads currently being conducted by the National 
Advisory Committee for Aeronautics with the relatively flexible swept- 
wing Boeing B-47 airplane have one phase of the program devoted to 
obtaining data where the controls are pulsed in a fashion to excite the 
aircraft structural freq^iencies. This report is concerned with presenting 
an analysis of some of the elevator-pulse data obtained in the Mach num- 
ber range from 0.62 to 0.82 as well as a reapplication of the results to 
more severe conditions . 

Since the frequency plane is well adapted to analysis of structural 
modes, it has been used in this report to determine the dynamic character- 
istics of the Boeing B-kj airplane and to demonstrate the effects of these 
characteristics on the pitching velocity and Incremental normal load 
factor. This type of analysis permits the simultaneous study of all the 
significant vibratory modes for the entire airplane. 
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airplane wing span^ ft 

airplane mean aerodynamic chords ft 

airplane lift coefficient^ L/qS 

rate of change of airplane lift coefficient with angle of 
attack per radian 

rate of change of airplane lift coefficient with elevator 
deflection per radian 
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rate of change of airplane lift coefficient with a per 
radian per second 

rate of change of airplane lift coefficient with pitching 
velocity per radian per second 

pitching -moment coefficient of airplane, M/q,Sc 

rate of change of airplane pitching-moment coefficient with 
angle of attack per radian 

rate of change of airplane pitching-moment coefficient with 
elevator deflection per radian 

rate of change of airplane pitching-moment coefficient with 
pitching velocity per radian per second 

rate of change of airplane pitching-moment coefficient 
with a per radian per second 

rate of change of airplane pitching-moment coefficient with 
lift coefficient 

acceleration due to gravity, 52.2 fty^sec^ 

pressure altitude, ft 

p 

airplane moment of inertia, slug-ft*^ 


index of to 

airplane radius of gyration about pitching axis, ft 


lift, lb 

airplane mass, W/g, slugs 
Mach number 

incremental airplane load factor 


dynamic pressure. 
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dynamic pressure in vicinity of horizontal tail, Ib/ft 
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St 

t 

T 

V 

Ve 

W 


wing area, sq ft 
h.orizont8Ll-tail area, sq ft 

time, sec 

period, sec 

true airspeed, ft /sec 


equivalent airspeed. 



airplane weight, Ih 

length from center of gravity of airplane to aerodynamic 
center of tail (negative for conventional airplanes), ft 


Ki, Kg,! 

K6 J 


coefficients of transfer function relating 0 and 6 g 
(see appendix A) 


a 

°t 


wing angle of attack, radians 
tail angle of attack, radians 


a 

^'f 

7 

9 

0 

Se 


rate of change of a with respect to time, radians/sec 
flight-path angle, radians 

dynamic response factor 
angle of pitch, a + 7 ^, radians 
pitching velocity, radians/sec 
elevator deflection, radians 


Tv = 


ot 

Sm5 


downwash angle, radians 
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r\^ tail efficiency factor^ 

p mass density of air, sliigs/cu ft 

D mass density of air at sea level, slugs/cu ft 

0 

I critical damping, percent 

D differential operator, d/dt 

s Laplace transform variable 

$ phase angle, deg or radians 

CO frequency, radians/sec 

parameter defined in gippendix A 


Subscripts : 


A 

E 

e 

eg 

F 

i 

n 

o 

E 

S 

t 


airplane 

elastic 

elevator 

center of gravity 

flexible 

input 

natural 

output 

rigid 

quasi-steady 

tail 

phase angle at co = 0 radians/sec 
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Matrix notation: 



square matrix 


column matrix 


For sign conventions used see figure 1. Dots over a symbol refer to 
derivatives with, respect to time, and a bar over a letter represents 

maximijm value. !Ehe absolute value of any term is denoted by I 1 . 


AIEPLAUE AUD INSTRUMENTS 


The airplane used in the tests Was a Boeing B-l^-TA aiiplane. A 
three-view drawing of the airplane is presented in figure 2 and perti- 
nent dimensions are given in table I. 

Because of the nature of the flight programs, it was necessary to 
install a large amount of instrumentation; however, only those instru- 
ments which were used to measure quantities pertinent to this paper will 
be mentioned in detail. 

Linear accelerations were measured by standard KACA air-damped 
recording accelerometers and by remote-sensing oil-damped Statham elec- 
trical accelerometers in conjunction with oscillograph-type recording 
instruments. Angular velocities and angular accelerations were measured 
about the three airplane reference axes by rate -gyro -type electrically 
differentiating, magnetically damped standard NACA tunmaeters. Control 
deflections were measured by remote -sensing electrical transmitters 
linked directly to the control siarface in conjunction with a galvanometer. 
The angle of attack was measured by a flow-direction recorder having a 
vane-type sensing device located on a boom 73 inches ahead of the nose. 

Electrical resistance-type wire strain gages were located through- 
out the structure. The output of these strain gages was recorded on 
oscillograph- type recorders. An NACA optigraph system for photographi- 
cally measuring deflections was also incorporated into the instrumentation. 

The data obtained from the instruments were recorded photographically. 
The recorders were synchronized at 0.1-second intervals by means of a 
common timing circ\iit. All instruments were damped to about O .67 of 
critical damping. The location of the pertinent instaruments or sensing 
elements is shown in figure 2. A simimary of the quantities measured. 


NACA IM L^hU-S 


Cl 


7 


instrument locations, types of Instrument, instmmient accviracies, and 
natural freqiiencies is given in table II. 


DATA-EEDUCTION ERRORS 


A detailed discussion of the factors affecting the accuracy of the 
Fourier method and frequency-response theory as applied to aircraft 
dynamic analyses is given in reference 3« This section, therefore, will 
be devoted to outlining briefly the main sources of error and the pre- 
cautions taken to minimize the errors affecting this paper. 

Two items which contribute the largest source of error in the data 
reduction are the reading accuracy and the reading time inteirval. In 
an effort to obtain maximum reading accuracy, time histories of the per- 
tinent parameters were read on the miversal telereader and telecordex 
in conjunction with an IBM card punch. With four-power magnification, 
the accuracy level of the reader is 'tO.003 inch of film deflection. The 
record sensitivities used and the accuracy levels (based on the reading 
acciiracy for the primary data coordinates) that they signify are as 
follows: 


Record sensitivity Accuracy level 


0 0.249 radians /sec /in. ± 0.00075 radian/sec 

A& l4 deg/in. t0.042 deg 

An 1.013g/in. ±0.005g 


These reading errors will result in an error in the frequency-response 
amplitude ratio proportional to 1 /cd which is substantial only in the 
low-frequency range. Records were read at a reading interval 
At = 0.02 second. In order to have a minimum of six data points defining 
each fxall cycle, a reading Interval of At = 0.02 second would pemlt 
a measured maximum frequency of 

CO = = _2jL_ = 52.4 radians/sec 

6 At 0.12 ' 

For a reading interval of At = 0.02 second, a rough estimate of the 
error in phase angle between the aircraft output and input for a given 
frequency (for example, co = 45 radians/sec) woiild be 

$ = CO Ab = 43 ( 0 . 02 ) = 0.86 radian = 49° 


This is a substantial error in the phase angle for the highest frequency 
band, and its effect is treated elsewhere in the paper by considering the 
total effect of phasing. 
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In transforming a ti^sient to the frequency plane by use of the 
Fotirier method, one of the important factors which must be considered 
is frequency content of the input. The pvilse method is desirable as a 
flight-test technique because it enables a frequency response to be com- 
puted from one flight run of a few seconds duration; however, it also 
has the limitation that the shape of the pulse input has an important 
effect on the harmonic content of the input. An input must have enough 
frequency content to excite svifficiently all harmonics over the consid- 
ered frequency range in order to provide an accurate frequency response. 
As the triangular pulse approaches a pure impulse, the frequency content 
of the triangular p\ilse approaches the ideal frequency content; however, 
as the length of the base increases, this frequency content will 

bottom (approach zero-frequency content) at a frequency of approx im ately 
0 ) = This condition produces an inacc'urate frequency response at 

the bottoming frequencies, (it should be noted that bottoming is a char- 
acteristic of symmetrical inputs.) Thus, frequency content of the input 
must be examined carefully in order that sufficient excitation be pro- 
vided over all the frequency range considered. As mentioned previously 
the pulse maneuvers used in this paper were selected in order to have 
good" frequency content over the frequency range investigated. 

Instrument errors are given in the section entitled "Airplane and 
Ihstruinents . " Corrections due to the dynamic response of Instruments 
are made in the frequency plane eifter the harmonic ana, lysis of the tran- 
sient has been made. 



In the usual open-loop diagram as shovn above an elevator input 6^^ 
goes through a control position recorder and is recorded as 8 q; the 
same control input 6^^ restilts in the aircraft response Uq which goes 
through an Instrument and is recorded as Iq. The apparent or measured 

frequency response without instrument corrections is then — (i^)) . Now 

®o 

if R 2 (io)) is the frequency-response data of the control position 
instrument and R 2 (lts) the frequency-response data of the response 
recording instrument, then the desired aircraft frequency response 

— (ioj) is determined from the relation 
^i 
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Rj(M 

R2(i<^) 


If the natioral frequencies of both of the recording instruments are close 
to each other, the amoimt of instrument correction to the frequency 
response is very small. An examination of table II indicates that this 
is the case for the instruments used on this aircraft. 


From the foregoing discussion it is apparent the accuracy level of 
the results is least in the very low frequency bands (less than 
0.75 radians/sec) in the very high frequency bands (above 40 radians/sec), 
and in the bands where the input frequency content tends to bottom. 


METHODS 


The term stinactural flexibility may have various interpretations 
depending upon the viewpoint. In this paper it refers to two separate 
phenomena, that of aircraft static aeroelastic deformations and that of 
deformations due to the structural or vibratory modes. References are 
frequently made and results given for the B-47A Boeing airplane with three 
separate amomts of flexibility; (l) an elastic or actual configuration 
with all its structural modes, ( 2 ) a quasi-steady configuration where 
response is affected only by quasi-steady aeroelastic effects such as 
would be measured in very slow maneuvers where higher vibratory modes are 
neglected, and ( 5 ) a rigid configuration whose characteristics are derived 
from the measured results. For the sake of brevity, these configurations 
are subsequently referred to as the elastic, quasi-steady, and rigid 
aijrplanes . 


Selection of Tests 

In this section, the methods used in selecting the flight runs for 
analysis and for reducing the results are outlined in a general way; addi- 
tional details are given in the appendix. 

The flight runs selected for detailed analyses were all short 
elevator-pulse maneuvers with the other controls held neutral. For these 
rims, the weight, altitude, and center- of -gravity position were nearly 
constant whereas the Mach number varied; however, some effects of altitude 
and center-of -gravity position on the aircraft short-period natural fre- 
quency were examined. Inasmuch as the effects of flexibility, transfer 
coefficients, and stability derivatives were all to be obtained from the 
data, use of the frequency plane was thought to be the most effective 
technique for evaluation of these quantities. As a result, only those 
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longitudinal maneuvers having elevator-pulse inputs of sxifficlently high 
frequency content to excite as many of the aircraft symmetrical struc- 
tural modes as possible were selected. Some typical time histories of 
the quantities used in the analyses of these maneuvers are given in 
figure 5 for three different Mach niimbers. 


Data Reduction 

In order to obtain the frequency responses of the airplane, the 
time histories, such as given in figure 5, were transformed to the fre- 
quency plane by means of the Foiirier method. The Foirrler integrals of 
inputs and outputs were evaluated by using the IBM equipment and the 
method of integration of product curves described in reference 4. The 
frequency responses were obtained by tedsing the ratio of the Fo\irier 
integraJ. of the output to the Fourier integral of the input as indicated 
by the expression 




An(t)e-i^‘fcdt 



A5e(t)e“^“'‘*^dt 


The frequency responses of the measured quantities 




f An 


felose’ W 


/eg 


'eg 


y and ( I thus obtained are shown in figure 4 


wing tip 


A6. 


for three different Mach numbers. 


In this analysis^ an attempt was made to extend the frequency range 
investigated as high as possible in order to get as many of the signifi- 
cant symmetrical vibratory modes as possible. In order to separate the 
short-period mode from the frequency- response curves of the pitching 
velocity^ points were computed at eveiy o) = 0.55 radian per second up 
to k radians per second. Over the rest of the frequency range ^ points 
were measured at frequency intervals of = 1.5 radians per second 

with additional fill-in points in the region of the structural modes. 

For the normal load-factor frequency-response curves the points were com- 
puted at frequency intervals of Acd = 1.5 radians per second with no 
fill-in points. In this case^ the intervals may not be sufficiently 
small to define the short-period modes as well as the measured pitching 
velocities. 

In figure 5 the predominant-symmetrical airplane frequencies are 
shown. The lines are the result of selecting the frequencies at which 
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the peak amplitude ratio in the Fourier analysis occurred. The test 
points were obtained from the results of a separate study of 100 runs 
involving I7OO traces in which the periods (in seconds) of all oscil- 
lations below 50 radians per second appearing on the recoils were con- 
verted to frequency and plotted against Mach number. The range of Mach 
number covered in this study was 0.14-5 ^ M < 0.814-; the range of center- 
of -gravity position was 0.12c to 0.50c; and the range of altitude was 
20,000 feet to 55^000 feet. 

In figure 6 the natural frequency of the short-period mode is shown 
plotted against Mach number for various altitudes and center-of -gravity 
positions. The lines are arbitrarily drawn fairing lines through the 
data points. 


Determination of Frequency Response of Rigid Conf lg\iratlon 

In order to extrapolate resiUts obtained from the flexible aircraft 
to those of the rigid aircraft, it was necessary to determine first the 
transfer coefficients and longitudinal stability derivatives of the flex- 
ible aircraft. The transfer coefficients in pitch were determined from 
the frequency-response data such as given in figure ^4- by the vector- 
least-squares method of reference 5* In this analysis, the short-period 
regime of the frequency-response curve was fitted by this vector- least- 
squares method and the transfer coefficients shown by the test points in 
figure 7 were obtained. The solid lines thro\jgh the data points for the 
plots of K]_, K2, and against Mach number are arbitrarily drawn 

fairing lines. Inasmuch as the Kg coefficient exhibited considerable 

scatter, the solid fairing line through the data was determined on the 
basis of the faired lines for the Kq_, K2, and coefficients and a 

value of Cr determined from flight measurements. These flight values 

of the lift-curve slope are shown in figure 8 and were determined 

from an analysis of 25 runs obtained during slow push-pull maneuvers by 
plotting airplane lift coefficient against angle of attack a. 

The stability derivatives shown in the remainder of figure 8 were 
obtained from the selected flight data and are those nonnally associated 
with the short period of the flexible aircraft which Includes static aero- 
elastic deformations. The computation of the derivatives of figure 8 was 
based on the equations of motion in the appendix, the valiies of A from 
reference 6, and the aircraft mass and geometric parameters listed in . 
table I. Also shown in figure 8 are the nat\xral frequency and percent 
critical damping which were obtained from the transfer coefficients and 
the equations shown in the appendix. 
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In order to evaluate the effects of aeroelastic deformations on the 
Boeing B-47A airplane, the stability derivatives and transfer coefficients 
were modified for the rigid configuration. In this modification, the 


flexible -to-rigid parameters 




shown in 


figure 9 were taken from reference 7 and based on tests of models 


having various degrees of stiffness. 


The parameter 




was 


taken frcan reference 8. These parameters of figure 9 are for a gross 
weight of 125,000 pounds and have negligible variation for the range of 
gross weight in this paper. These parameters were determined on the basis 
of changes in the spanwise lift distribution, a shift in the aerodynamic 
center, changes in the tall angle of incidence, changes in control effec- 
tiveness, and so forth which were caused by static deformations of the 
aircraft struct\ore due to flexibility in a swept-wing configuration. The 
parameters of figure 9 iu conjunction with the relations contained in the 
appendix were used to compute the stability derivatives and transfer coef- 
ficients of the rigid airplane. The stability derivatives of the rigid 
conflgiiratlon and transfer coefficients thus computed are Indicated in 
figures 7 aud 8 by broken lines. 


The frequency response in pitching velocity for the rigid airplane 
shown in figure 10 was then computed by using the rigid transfer coef- 
ficients and frequency-response relations given in the appendix. Also 
shown in figure 10 are the frequency responses in pitching velocity for 
the elastic airplane and quasi-steady airplane. The frequency response 
of the elastic airplane was shown previously in figure 4 and was computed 
directly from the measTired flight data. The frequency response of the 
quasi-steady airplane is the short-period mode of the elastic aii^Jlane 
with the structural modes removed. 


D7MMIC-RESP0NSE AMLYSIS 


Most of the figures presented thus far have dealt with data and 
parameters extracted from the data as obtained from the fli^t-test 
maneuvers. In order to study more fully the effects of flexibility on 
the loads and longitudinal responses of the aircraft, however, further 
computation with these parameters and the extension of the data to more 
hazardous conditions are required. !I3ie response characteristics of the 
elastic aircraft used in combination with a severe but practical type of 
control input permit extension of the data beyond the test range. 

The triangular pulse method used in reference 9 in a study of the 
mn mum dynamic loads during landing appeared to be the most appropriate 



MCA RM L54L16 


15 


type of Input and was selected on the hasis of several desirable fea- 
tures: The frequency content of the triangles could easily be varied by 

changing the time base of the triangle; the triangles had sufficient 

frequency content to excite the higher airplane vibratory modes; Eind the 
triangular input was well adapted to synthesizer operation and loads com- 
putations. The triangles used in this investigation were Isosceles. 

In some of the computations of response to triangular Inputs, the 
electromechanical Fourier synthesizer was used to produce transient- 
response curves from frequency- response data. In other cases the com- 
putation was carried out on calculating machines. The operation of the 
Fourier synthesizer is described in some detail in reference 

Figure 11 shows three computed time histories of the response in 
pitching velocity of the elastic, quasi-steady, and rigid airplanes to a 
triangular input. In this case, the results are shown for a period ratio 
of Tj^^T = 0.1 and a Mach number of 0.82. In the computations of fig- 
ures 11 to 20, the maximum elevator deflection is 0,1 radian. 

By varying the ratio Tj_/T of the input over a range of values, a 

series of plots similar to those of figure 11 were obtained. The maximum 
value of pitching velocity of the elastic airplane following the trlan- 

• # 
gular piilse was designated as 0^^ of the quasi-steady airplane as 

and of the rigid airplane as 0p. A comparison of these maximum transient 

responses was then made by plotting their ratios against the period 

ratio Tj[^T as shown in figure 12. The ratio Indicates the 

effect of static aeroelastic deformation, the ratio 0 e/ 9s indicates the 

effect of the structirral modes, and the ratio ^e/®r indicates the over- 
all effects of flexibility on the longitudinal response in pitching veloc- 
ity of the airplane. 

The data of figure 13 are a conversion of the ratios of figure 12 to 
a percent alleviation or magnification factor of the pitching velocity. 

The ordinate is constructed so that positive values indicate alleviation 
and negative values indicate magnification. 

Plots similar to those of figure 11 were also obtained for the load 
factor at the center of gravity. A plot of the ratios of maximum load 
factor similar to those of maxim\am pitching velocity of figure 12 was 
obtained and is presented in figure l4. The same trend in normal load 
factor as in pitching velocity is indicated. As noted previously, the 
ratio Ang^An-p indicates the effects of static aeroelastic deformation; 
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the ratio /iaa^Aiig indicates the effect of structural modes; and the 
ratio Aug/zhip indicates the overall effects of flexibility on the nor- 


mal load-factor response. 

Another of the effects of flexibility obtained from figure 11 is 
the change in the time to reach peak pitching velocity for the elastic, 
quasi-steady, and rigid airplanes. Figure 15 is a plot of the time to 
reach peak pitching velocity against the period ratio for the 

three alip)lanes. 

A similar plot of time to reach peak load factor at the center of 
gravity against the period ratio for the elastic, quasi-steady, 

and rigid airplane is shown in figure l6. 

In the study of flexibility it is of interest to assess the impor- 
tance of the phasing relation between the modes. Up to this point, all 
results have included phasing but, in order to study the importance of 
this phase relation, the calculations were repeated with phasing removed 
and the results compared with the previous computations. Phasing was 
removed from the computation by considering the phase angle of the air- 
plane as constant and equal to the value of o) = 0 over the entire fre- 
quency range. The ratio of maximum pitching velocity, not including 
phasing, to the maximum pitching velocity, including phasing, plotted 
against the period ratio is shown in figure 1?; whereas, figure l8 is a 
similar plot for nmylTmini incremental normal load factor at the center of 
gravity. Inasmuch as considerable errors can occur in the phase-angle 
computation in the high-frequency range of the spectrum, these compari- 
sons also aid in determining the effects of such errors on the maximxmi 
loads and longitudinal responses. 

Another important parameter in the study of loads associated with 
flexible stiTictures is the dynamic -response factor (see ref, 9) defined 

^^^dvn 

as 7 = ^ — where ^sn^yn refers to the maximum amplitude of Incre- 

static 

mental load factor and ^static refers to the steady-state amplitude 

of the Incremental load factor for the same input. The variation of the 
dynamic response factor j is shown in figure 19 plotted against the 
period ratio T^^T for the elastic, quasi-steady, and rigid airplanes 

at a Mach number of 0.82 and for the elastic and quasi-steady airplanes 
at a Mach number of 0.66. 


In order to obtain an overall pictvure of the deformations and vibra- 
tions taking place in a typical flight maneuver, a three-dimensional plot 
of the time histories of the wing and fuselage bending defoimations d\aring 
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a pulse maneuver is shown in figure 20. This is the same maneuver as 
shown in figures 3(a) and ^(a). 


DISCUSSION 


The methods employed in this paper are all based on the assumptions 
of linearity and the small perturbation theory. The flight-test tech- 
niq,ue of pijlsing the elevator is essentially a small-perturbation method. 
Inasmuch as the flight conditions and maneuvers analyzed were selected 
so that the airplane would respond in the linear range of the transfer 
coefficients and stability derivatives^ the qiiantities obtained from 
this analysis apply most accurately to this linear range. 

The frequency- response curves shown in figure k indicate the dynamic 
characteristics of the elastic airplane such as the short-period mode and 
the vibratory or structural modes which appear as peaks on the amplitude- 
ratio cu 2 Tves. Evidences of the airplane phugoid mode do not appear as 
is to be expected because the phugoid is a long-period mode which cannot 
be detected from the data with pulse-type inputs. In general^ a measure 
of the relative importance of each of the individual modes on the response 
of the airplane is given by the resonant amplitude of the modes and by 
the frequency at which the modes occur. The higher the frequency at which 
the mode occurs, the less effect the mode will have in the normal airplane 
maneuver range since any practical pilot-applied forcing function will 
excite higher harmonics with lower amplitude. The structural mode at 
03 = 27-5 radians per second appears as a dominant mode in all the 
frequency- response traces; this is also particularly noticeable in the 
three-dimensional plot of airplane deflections shown in figure 20. 

The plot of symmetrical airplane frequencies shown in figure 5 indi- 
cates an increase in the natural frequency of the short-period mode with 
increasing Mach number and, as expected, little or no change with Mach 
number of the structural mode frequencies. Mach nimiber, however, has a 
definite effect on resonant amplitudes of the structural modes and gen- 
erally shows an increase in amplitude with increasing Mach number. An 
attempt was made to identify the structural frequencies in figure 5 vith 
the usual uncoupled cantilevered modes. This attempt proved imsuccessful 
inasmuch as the airplane vibratory modes are coupled modes containing 
elastic and aerodynamic coupling and, because this is a swept-wing con- 
figuration, each wing bending mode has its own associated torsion. 

In figure 6 the effects of either a rearward movement of the center- 
of-gravity position or an increase in altitude is seen to decrease the 
natural frequency of the short-period mode in accordance with expectations. 
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The effects of static aeroelastic deformations in the airplane 
transfer coefficients and stability derivatives are indicated in fig- 
ures 7 and 8 hy the difference between the solid and broken lines. The 
force derivatives and transfer coefficients for the flexible air- 


a 

plane are shown to be less than those for the rigid airplane. The ele- 
vator effectiveness derivative CL is reduced because of flexibility. 

iiig 


The largest effect due to flexibility is evidenced in the derivative 


From the curves of short -period natirral frequency and percent 

critical damping | shown in figure 8 for the flexible and rigid air- 
planes it is evident that flexibility reduces the short-period natural 
frequency and increases | . The reduction in appears to be due to 

the reduction in CL • 


The frequency responses shown in figure 10 for the three different 
amounts of flexibility in the airplane show that^ as elasticity is intro- 
duced into the rigid structure^ a change takes place in the frequency 
response as evidenced not oixly by the change in the short-period mode but 
also by the appearance of vibratory modes. 


In the typical time responses of the elastic, quasi-steady, and rigid 
airplane shown in figure 11, the main change in the time response of an 
airplane for three separate amounts of flexibility is seen to cause the 
elastic-airplane response to oscillate around the quasi- steady-airplane 
response and to modify the maximum amplitudes of the responses and times 
to peak amplitude for each case. 


Some pulses applied to the airplane Imposed more severe rates of 
elevator displacement than could be achieved by a pilot in flight. An 
analysis of the pulse maneuvers of the flight-test program indicated a 
maximum elevator rate of 0.75 radian per second which is equivalent to 
Tj^ 0.27 second for a 0.1- radian deflection. For the calc-ulations pre- 
sented in this paper, this value represents a period ratio ^± 1 ^ = 0.1 

which is the minimum that the pilot can apply to the airplane by control 
inputs. Under unusual conditions, however, when the forcing function is 
attributable to gusts or separation phenomena, smaller period ratios may 
be attained. In this lower range, however, the computed values may be 
high because unsteady lift effects were not considered in the calculations 


The ratio figure 12 indicates that static aeroelastic 

deformation in all cases reduces the maximum airplane response in pitching 


velocity. The ratio Qe/qs indicates an increase in maximum pitching 
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velocity over the low ranges of Tj_^T where the structiaral modes are 
excited and indicates no effect over the high ranges of where 

there is insufficient frequency content in the input to excite the struc- 
tural modes • The ratio ®e/®e indicates the overall effects of flex- 
ibility on pitching- velocity response. Because of the excitation of the 
structural modes in the lower range of the maximum response in 

pitching velocity is increased. In the higher range of where the 

structural modes are not excited, static aeroelastic deformation decreases 
the maximum airplane response in pitching velocity. 

The results given in figure 15 show that static aeroelastic defor- 
mation provides approximately 11 -percent alleviation over the rigid case, 
and the stinctural modes provide magnification of the response up to 
90 percent. 

The ratios of maximum incremental load-factor values at the center 
of gravity shown in figure Ih indicate the same trends as the maximimi 
pitching- velocity responses, that is, an alleviation due to static aero- 
elastic deformation and a magnification due to the structural modes when 
excited. 

The time to reach peak pitching velocity as shown in figure I5 is 
a minimum for the rigid airplane and a maximum for the quasi-steady air- 
plane; the time for the elastic airplane varies between the two limits. 

The effect of static aeroelastic deformation is therefore one of 
increasing the time to reach peak pitching velocity over the rigid case; 
whereas, the effect of the structural modes is to decrease the time to 
reach peak pitching velocity over the quasi-steady airplane. 

Ihe effects of the structiiral frequencies on the time-to-peak-load 
factor is shown in figure 16 to be a decrease in the time-to-peak-load 
factor when the structural frequencies are excited. 

The effect of removing the phase-angle relation in the computations 
as shown in figure I7 is to increase the maximimi pitching-velocity 
response in all cases. This is in accordance with expectations inasmuch 
as, in a simple vibratory system, the effect of removing the phase-angle 
relation is equivalent to removing damping from the system. 

The effect of excluding the phase-angle relation from the computation 
of maximum load factor at the center of gravity as shown in figure I8 
indicates the same trend as for the case of maximum pitching- velocity 
response, namely, to increase the maxlm’um load factor in all cases. 

The dynamic response factor 7^ of the quasi-steady airplane is 
approximately 1.2 for a Mach number of 0.82 and I.3 for a Mach nmber 
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of 0.66. (See fig. 19* ) The curve for dynamic response factor y„ 

* E 

against the period ratio for the elastic airplane exhibits two 

definite peaks, one due to the structural modes and the other due to the 
short-period mode. The maximum dynamic overshoot of the elastic airplane, 
however, is the same as the quasi- steady-airplane dynamic overshoot. It 
should be noted that, if the short-period natural frequency and vibratory- 
mode frequencies were closer together, then the dynamic response of the 
short-period mode and vibratory modes could add and form a critical 
loading condition. (The dynamic response factor 7^^ for the rigid air- 
plane is approximately 1.43 for a Mach number of 0.82.) The overall Mach 
number effect appears to be a decrease in dynamic response factor 7 
with increasing Mach number due in part to an increase in ^static 'wi'^h 
Mach number. 

The three-dimensional time histories of wing and fuselage defor- 
mations shown in figure 20 indicate that all the modes are coupled. The 
fuselage deformations can be seen to be the superposition of the vibra- 
tory modes on the static aeroelastic deformation. 


CONCLUSIONS 


Effects of structural flexibility have been examined by considering 
the swept-wing Boeing B-47A airplane with three different amomts of 
flexibility and comparing the maximum responses of each to the same inputs. 
The analysis of some of the elevator-pulse data of this paper apply spec- 
ifically in the Mach nimber range from 0.62 to 0.82 in the linear range 
of longitudinal response, and the following conclusions are established: 

1. Static aeroelastic deformation alleviates loads and decreases 
the maximum pitching- velocity response. 


2. The structural modes ■when excited increase loads and maxim-urn 
pitching- velocity response . 


3. The magnitude of the force derivative C_ , two moment deriva- 

•^a 

tives C_ and CL , and the transfer coefficients are less for the 
"quasi -steady" airplane than for the "rigid" airplane. 


4. The short-period natural frequency of the "elastic," "qxaasi- 
steady, " and "rigid" airplane decreases with increasing altitude, with 
rearward movement of the center of gravity, and with decreasing Mach 
number. 
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5 . The structural frequencies of the "elastic" airplane shov little 
change with Mach number or altitude. 

6. The "quasi-steady" airplane has a lower natural short-period fre- 
quency 0 ^ Sind a somewhat higher percent critical damping | than the 

"rigid" airplane. 

7 . The inclusion of phasing in the computation of longitudinal coor- 
dinate responses and loads of the airplane reduces the maximum values. 

8. The time-to-peak-pitchlng-velocity response and Incremental load- 
factor response of the "quasi -steady" airplane is higher than that for 
the "rigid" airplane or the "elastic" airplane when the structural modes 
are excited. 

9 . Airplane structural frequencies should be well separated from 
the short-period mode and should be sufficiently high to be beyond the 
pilot's ability to excite them with control inputs. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., December 8, 195^. 


20 


MCA EM L^hLLG 


APPENDIX 

SUMMARY OF MA03IEMATICAL RELATIONS USED IN EVALUATING DATA 


The longitudinal eqiiations of motion for horizontal flight used in 
this analysis (see ref. 2) are 

X to = Ci^ + CLj /iBe (Al) 


^ e = to. + + c„.e + (A2) 

and 

9 =a+7^=a+^An (AJ) 

These equations are "based on the usual assumptions of linearity, small 
angles, and no loss in airspeed during the maneuver. 

The usual transfer function between pitching velocity and incremental 
elevator deflection o"btained by solving equations (Al) to (A5) simulta- 
neously is given "by 


D 0 + K]_D9 + K20 = K^DSg + KgSg 


(a4) 


where the transfer coefficients K]_, K2, K^, and Kg are defined as 

follows : 


K _ <lS 

^1 - Y~ 




Ct 


Kg = - + Cme 


(A5) 

(a6) 


K5 = ^ CLgC,. 


(a7) 
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f Ct 

I mV \ “6 



(A8) 


The pitching-velocity transfer coefficients K^, nnii 

of equation (a 4) were computed from the pitching -velocity frequency- 
response data hy use of the vector-least-squares method of reference 5 
as follows ; 



and N is the number of data points used. 


The following expressions and procedures were used in the deter- 
mination of stability derivatives. G?he Ct derivative, as mentioned 

■“(X 

in the text, was determined from a separate investigation of slow push- 
p\iU maneuvers. Values of the ratio A used in this report were 
obtained from reference 9* The remaining longitudinal derivatives were 
computed by use of the following approximate relations from reference 2: 
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+ Cni.'] 

aj 

1! 

1 

9 

(A12) 



(A13) 


"qp" ^ ^ *^0 

(A1J+) 


I ^5 

(A15) 

<lSc 1 _ qg c_ 

mV ^ ^0 


^Lg = 5E^ ^5 

(A16) 


tl 

(A17) 


Kt 

i = — 

2o^ 

(A18) 


In converting the rigid stability derivatives from tbe derivatives 
of the flexible airplane, the following definitions and expressions were 
used: 


(%), ' 

(A19) 

II 

(A20) 



Wr /acd \B6 /y 

(A21) 

\^/f 
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(A22) 


(A25) 


Frequency response for the rigid airplane in pitching velocity was 
coDiputed hy using the following relations for amplitude ratio and phase 
angle; 



(A2l<-) 


(A25) 
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TABLE I.- PERTmENT PHYSICAL CHARACTERISTICS AMD 

DIMEMSIONS OF BOEING B-^TA TEST AIRPLANE 

Total wing area, sq ft 

Wing spaji, ft U6 

Wing aspect ratio 9 >^3 

Wing thickness ratio, percent 12 

Wing taper ratio 0.42 

Wing mean aerodynamic chord, in 155-9 

Wing sweephack (25-percent chord line), deg 55 

Total horizontal-tall area, sq ft . 268 

Horizontal-tail span, ft 55 

Horizontal -tail mean aerodynamic chord, in 102.9 

Horizontal- tail sweephack (25-percent chord line), deg 55 

Distance from center of gravity to quarter chord of horizontal- 

tail MAC (center of gravity at 0.21c), ft 4? 

Airplane weight, lb 125,000 to 127,000 

Airplane pitching radius of gyration, ky, ft 19*5 

A 0.18 



TABLE II.- SUMMARY OF INSTRUMENTATION AND ACCURACIES 


ro 

G^ 


Quantity measured 

Where measured 

Instrument 

used 

Units 

Instrument 

range 

Instrument 

accuracy 

Natural 

frequency, 

radians/sec 


Center of 
gravity 

Three- component 
air-damped 
accelerometer 

g 

-1 to if 

0.0125 

66 

Nonnal accelera- 
tion 

Tail 


g 

-2 to 6 

0.0200 

66 

Optlgraph 

housing 

Two- component 
air-damped 
accelerometer 

g 

-1 to 4 

0.0125 

66 

10 viog stations 

Statham electrical 
accelercmeter 
and galvanometer 

g 

il2 at tin to 5*05 at inhoard station 

^ -4.25 

Recording 0.1 

Accelercmeter 628 
Recorder 75 

6 fuselage 
stations 

do 

g 

to 5.8 
-2.51 -7.8 


Do. 

Pitching 

velocity 

Center of 
gravity 

Angular velocity 
recorder 

radians 

sec 

±0.25 

0.005 

85 

Pitching 

acceleration 

Center of 
gravity 

Angular accelera- 
tion recorder 

radians 

flec^ 

±0.50 

0.010 

44 

Elevator 

angle 

Midelevator 
semi span 

Control position 
transmitter and 
galvanometer 

deg 

25 up and 15 down 

0.40 

75 

Angle of 
attach 

73 inches in 
front of nose 

Vane -type trans- 
mitter of flow- 
direction 
recorder 

deg 

±30 

0.10 

94 

Periods of vari- 
ous structural 
frequencies 

Several vlng 
and fuselage 
locations 

Strain gages and 
accelerometers 

sec 

Varied 


75 

Time 


Timer 

sec 


Approximately 0.005 


1 




wa vovM 
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Figure 1.- Sign convention employed with positive directions shown. 



116* at 0°dihedral- 
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o standard accelerometer 
□ Statham accelerometer 
O Angular velocity and angular 
acceleration recorder 
A Control position transmitter 

* Strain-gage bridge 

• Optigraph torget 





Optigraph housing 


Figure 2.- Three-view drawing of test airplane, 



Pitching 

velocity. Incremental 

Incremental load factor, An, g units A radians elevator angle, 


NACA RM L5ALI6 



(a) M = 0 . 82 . 

Figure 5*- Time histories of measured quantities. Hp = 50,000 feet 
center-of- gravity position at 0 . 21 c. 


CM 



Pitching 

velocity, incremental 

Incremental load factor, An, g units a radioni elevator angle, 

’ » second a 6. deg 
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Time,t, sec 


(b) M = 0.75. 
Figure 5»- Continued. 
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(c) M = 0.66. 


Figure 5*- Concluded. 
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Mach number, M 


Figure 5.- Predominant symmetrical airplane frequencies. 
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Frequency, CO, Frequency, CO, 


36 


MCA EM L54.L16 



Figure 6.- Effect of altitude and center-of-gravity position on short- 
period natural frequency of test airplane. 
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Figure 7 .- Variation of pitching-velocity transfer coefficient with Mach 
mimber for test airplane. Ep = 30,CX)0 feet; center-of- gravity posi- 
tion at 0 . 21 c. 
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airplane. Hp = 30,000 feet; center-of-gravity position at 0.21c. 
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Mach number, M 


(Td) Moment derivatives. 
Figure 8.- Continued. 
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200 300 400 

Equivalent airspeed,^, 


Figure 9*“ Parameters used for transforming derivatives of elastic air- 
plane to those of rigid airplane obtained from references J and 8 for 
gross weight of 125,000 pounds. 
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Figure 10.- Comparison of typical elastic, quasi-steady, and rigid 
pitching-velocity frequency responses. ^ = 30,000 feet, center- 
of-gravity position at 0.21c, 
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Figure 11.- Typical pitching-velocity time-history responses to frequency 
responses of elastic, quasi-steady, and rigid configurations due to 
triangular-elevator-pulse input . Tj[/r = 0.1; M = 0.82; Hp = 30,000 feet 

center -of -gravity position at 0.21c. 
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Figure 12.- Variation of ratios of maximum pitching- velocity responses with 

period ratio. 
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(a)lncremental load factor ratios at M=.82. 
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(b) Incremental load-factor ratio, =■ ,at 
two Mach numbers. ^*^8 


Figure 1^+.- Variation of ratios of maximum incremental normal load factor 

with period ratio. 
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(a)M=82. 


-o Elastic 

□ Quasi-steady 


□ Quasi- steady 



Period ratio,- 


Figure l8.- Variation of ratio of maximum incremental normal load factor 
not including phasing maximum incremental normal load factor 

including phasing An with period ratio. 










(a) M = 0.82 and An = I.50. 

static 



(b) M = 0.66 and Ang^o^^^g = 0.64. 


Figure 19«- Variation of dynamic-response factor 7 with period ratio 
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